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Abstract. The sensitivity of future gravitational wave (GW) observatories will 

^.^ be limited by thermal noise in a wide frequency band. To reduce thermal noise, 

ly^ the European GW observatory Einstein GW Telescope (ET) is suggested to use 

(^ crystalline silicon test masses at cryogenic temperature and a laser wavelength 

* r^ of 1550 nm. Here, we report a measurement of the optical loss in a prototype 

O , high-resistivity crystalline silicon test mass as a function of optical intensity at 

Q room temperature. The total loss from both the bulk crystal and the surfaces 

• was determined in a joint measurement. The characterization window ranged 

^ from small intensities below IW/cm^, as planned to be used in ET, up to 

^ 21kW/cm^. A non-linear absorption was observed for intensities above a few 

^ kW/cm-^. In addition we have observed an intensity-independent offset that 

^^ possibly arises from absorption in the crystal surfaces. This absorption was 

<*^ estimated to Ogurf ~ SOOppm/surface, which might be too high for a cryogenic 

O , operation of a fibre-suspended silicon test mass. Such an offset was not observed 

I I in other recent measurements that were insensitive to surface absorption. Finally, 

a set of further characterization measurements is proposed to clearly separate the 

^— H contributions from the surfaces and the bulk crystal. 
> 

. 1. Introduction 

■^ At present the construction phase of the so-called Advanced GW detectors [l}[3] 

j^ nears its completion. The establishment of a gravitational wave (GW) astronomy, 

. however, requires a sensitivity improvement by another order of magnitude over 

• • the entire observation band. This motivated a design study for an European 3rd- 

^J^ generation GW observatory, the Einstein GW Telescope (ET) l4l, as well as first 

K/ concrete considerations for upgrades of the Advanced LIGO detectors [s] . The main 

^ sensitivity limit of the advanced detector generation at mid-frequencies will be thermal 

noise that originates from thermally excited motions in the mirror substrates, the 
dielectric reflective coatings and the suspensions. To overcome these noises, it was 
proposed to cool the test masses of the ET low-frequency interferometer to cryogenic 
temperatures and to employ crystalline silicon as test mass material [4] . The advantage 
of crystalline silicon is that it is available in the dimensions required for the ET test 
masses. A "drawback" is that with an absorption coefficient of a « 10 /cm at 1064 nm 
it is completely opaque at the wavelength that is currently used in GW detectors. 
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The employment of crystalline silicon will thus enforce a change of the operation 
wavelength, for example to 1550 nm |4 that is located within the silicon absorption 
gap. 

While predictions basing on photocurrent measurements led to the hope for sub- 
ppm/cm absorption at 1550 nm (6Jl7], to the best of our knowledge only recently first 
direct sensitive measurements of the optical absorption of undoped samples at 1550 nm 
were presented [8|9]. In [8] we have reported an optical absorption of a = 264ppm/cm. 
This value is two orders of magnitude higher than expected from the residual doping of 
the silicon crystal (the specific resistivity was specified to 11 kSlcm). This measurement 
comprised both the absorption in the bulk material as well as the surface absorption. 
In contrast, measurements done by the LMA group in France on a ~ 30kJ7cm crystal 
that used a beam deflection method ^ have found the absorption value of just the 
bulk crystal to be much lower in the order of < 5ppm/cm [9]. Furthermore, the 
measurements from the LMA group indicated that the optical absorption increases at 
high optical intensities. 

Here, we report a measurement of the optical loss in a prototype crystalline 
silicon test mass as a function of optical intensity. The measurement was performed 
in a silicon crystal that was highly reflectively (HR) coated at its surfaces and thus 
formed a monolithic cavity (Fig. [l](a)). Thus, the absorption in the bulk crystal, the 
dielectric coatings as well as in the surfaces beneath the coatings could be measured 
simultaneously as illustrated in Fig. Il](b). The characterization window ranged over 
more than 4 orders of magnitude from small intensities below 1 W/cm^ , as planned 
to be used in ET, up to 21kW/cm^. Our measurements confirmed the occurrence of 
intensity-dependent effects for intensities above a few kW/cm^. 

The high intensity-independent offset of ~ 250 ppm/cm observed in the 
measurements reported here can't be explained by the absorption in the bulk crystal. 
For a sample with a specific resistivity of llkficm this free-carrier absorption is in 
the order of a few ppm/cm as shown in Fig. p] This order of magnitude has also 
been observed by the LMA group, where a bulk absorption of < 5 ppm/cm was 
measured [9]. Furthermore, the contribution of the dielectric Si02/Ta2 05 coatings 
to the total absorption is with about 25 ppm/coating in a first-order consideration 
negligible fll] . This implies that the optical absorption of crystalline silicon needs to 
be treated separately as (a) the absorption in the bulk crystal and (b) a much higher 
absorption of about 800 ppm per transmission in the crystal surfaces. 

2. Method summary 

We performed absorption measurements on a silicon crystal at a wavelength of 
1550 nm. First absorption results measured with an optical intensity (here and 
in the following intensity refers to the peak intensity of the Gaussian laser beam 
power distribution inside the cavity) of 2.2kW/cm^ were published in ['8 . The 
cylindrical silicon crystal (shown in Fig. [l](a)) that had a diameter of 100 mm 
and a length of 65 mm was manufactured with the Czochralski technique. The 
specific resistivity was specified by the manufacturer to about llkJ7cm and was 
due to a residual boron contamination, thus meaning an impurity concentration of 
approximately 2 x 10^^ atoms/cm'^ 12]. The crystal end surfaces were polished to 
be convex curved with a radius of curvature (ROC) of 1 m and coated with a highly 
reflective Si02/Ta205 coating using ion beam sputtering. The power reflectivity at 
the wavelength of 1550 nm was measured to Ri = R2 = 99.96. Thus, a monolithic 
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Figure 1: (a) Photograph of the prototype test mass, the diameter is 100 mm. (b) 
Schematic of the monohthic cavity. Regions where optical absorption takes place, 
namely the bulk crystal, the polished crystal surfaces and the HR dielectric coatings, 
are highlighted. Independent measurements have shown that the contribution of the 
HR coating to the total absorption is, in a first-order approximation, negligible n][l3 



(c) Laser power reflected by the cavity (orange dotted line) and fit of the resonance 
peak (solid blue line) for an intensity of 0.4 W/cm^. Without the optical loss the cavity 
would be impedance-matched and the minimum of the resonance peak would drop to 
zero, since the HR-coatings have the same reflectivity (confirmed in an independent 
measurement). 



optical cavity with a free spectral range (FSR) of 663 MHz and a finesse of w 1500 
was available for absorption measurements. 

In 8 1 we have shown that the optical absorption per round-trip and the total 
cavity round-trip loss (Lrt) were identical within a few percent, which corresponded 
to the measurement error bars. Thus, a simple measurement of the light's impedance 
matching to the cavity (with known mirror refiectivities) provides the optical 
absorption per round trip. For this measurement, a 1550 nm laser beam was matched 
to the eigenmode of the monolithic cavity and detected in refiection. The laser 
frequency was linearly varied in the vicinity of a cavity resonance. Frequency markers 
imprinted with an electro-optic modulator provided an accurate calibration of the time 
axis. A more detailed description of the measurement setup is given in [8] . 

The laser power injected into the cavity was varied from 5 |i.W to 700 mW, being 
equivalent to intra-cavity peak intensities of 0.4W/cm^ -21kW/cm^, respectively. 
For each input power we performed 10 measurements of the cavity resonance peak 
using different scanning frequencies to ensure that the measurement was not limited 
by time-dependent effects. The peaks were fitted using the software developed for 
the photothermal self-modulation technique [14] and the impedance mismatch was 
determined (please note that the photothermal self-modulation technique itself was 
not applied for the measurements reported in this paper since at optical powers < 
1 kW/cm^ no absorption-induced deformation of the resonance peaks that is required 
for the method took place). The fitting parameters were the cavity in-coupling power 
reflectivity Ri and the effective out-coupling reflectivity i?2 = R2 + Lrt- R2 accounts 
for the entire cavity round-trip loss except for the transmission of the in-coupling 
coating. From this fit the round-trip loss was derived, the results are shown in Fig. [2] 
and discussed in the following section. Further geometric and material parameters used 
in the simulation were the crystal length of 65 mm and the silicon index of refraction 
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Figure 2: Optical absorption in our test mass derived from the impedance mismatch 
measurement as a function of the intra-cavity peak intensity. Left y-axis: absorption 
values scaled to the round-trip length. Right y-axis: total cavity round-trip loss. 



nsi = 3.48 [15| . Figure ^(c) shows an exemplary measurement at an intensity of 
0.4W/cm^ and a fit of the resonance peak (orange dotted line and solid blue line, 
respectively) . 



3. Discussion 



Figure [2] shows the results obtained for the optical absorption and a least-squares 
fit {y = 3.2 • 10~^a;^ -|- 2.7 • 10~^a; -I- 251.4) of the measurement distribution (orange 
dots and solid blue line, respectively). The left y-axis provides absorption values 
scaled to the round-trip length of 13 cm (thus in ppm/cm), while the right y-axis 
gives the total round-trip loss. Following the discussion given in [S], a maximal error 
of 10%, resulting from uncertainties in the simulation input parameters, the cavity 
geometric parameters and the mode-matching was assumed. As a consistency check, 
the measurements reported here were compared to the absorption value obtained for an 
intensity of 2.2kW/cm^ using the photo-thermal self-phase modulation technique [s] 
(shown in trace (c) of Fig. l3|. Both results excellently agree within the error bars, 
giving each time a value of a = 264 ppm/cm. 

The results of the measurement shown in Fig. [2] are twofold. First, our 
measurements show the impedance matching of the cavity to change in a non-linear 
way as a function of optical intensity. Assuming that this effect is due to bulk 
absorption, our results thus qualitatively confirm the measurements reported in [9] 



by the LMA group that used a beam deflection approach 10 , although the non- 
linear increase of the absorption was observed at an intensity lower by more than one 
order of magnitude in the LMA measurements. Degallaix et al. explained the non- 
linear behaviour by two-photon absorption processes. Such two-photon absorption 
creates free carriers that strongly increase the loss due to the free-carrier absorption 
when compared to just the absorption due to the residual dopant concentration. The 
free-carrier density and thus the effective absorption, however, depends on the free 
carrier lifetime and hence on the silicon crystal purity. As Degallaix et al. have shown 
in [9] , samples with higher purity (and thus a longer free-carrier lifetime) demonstrated 
a higher absorption for the same intensity. While the LMA measurements used a 
float-zone purified crystal with a specific resistivity of « SOkl^cm, our crystal was 
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Czochralski-grown with a specific resistivity of about llkJ7cni. Altliougli a detailed 
analysis of the free-carrier absorption as a function of specific resistivity is required, 
it seems likely that the disagreement in power dependence is due to the difference in 
the free carrier lifetime. 

Second, the absorption results of > 250ppm/cm reported here and in Js] are 
significantly higher than the results from other groups as shown in Fig. l3[(a,d). 
The measurement approaches employed there consisted in (i) spectral response 
measurements on solar cells [6] and (ii) beam deflection measurements fol. The 
first case only provides an indirect measurement, its results are therefore not 
unconditionally transferable to optical absorption and may be rather used to derive a 
lower limit as discussed in [s]. While the beam deflection approach is generally much 
better suited to characterize the optical absorption, it only provides information about 
the volume where the probe beam interacts with the thermal lens created by a strong 
pump beam. In [o] this region was selected to be within the optical substrate. Thus, 
merely the bulk absorption was evaluated, while no full information about the optical 
absorption of a complete test mass including both the bulk crystal as well as the 
surfaces was obtained. 

In contrast, our measurement approach optimally reflects the real application in a 
GW observatory, because it provides the full information about the optical round-trip 
loss. Both the bulk absorption and the absorption in the surfaces beneath the HR 
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Figure 3: Comparison of absorption values of crystalline silicon in the near infrared 
band, (a) Band-band absorption derived from photocurrent measurements [6J^. (b) 
Free-carrier absorption calculated for the residual doping in our sample (p-doping, 
N = 2 X 10^^ /cm'^) using the theory presented in 16 . (c, red circle) Total optical 



absorption in a prototype crystalline silicon test mass measured at an intensity of 
2.2kW/cm^ and scaled to the cavity round-trip length M. (d, head-down blue 
triangle) Bulk absorption of a silicon crystal measured with a beam deflection approach 
at an intensity of 2.2kW/cm^ 19]. 
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coatings contribute to this loss as illustrated in Fig|I](b). When compared to the LMA 
measurements, it stands out that at low intensities our data converges to an absorption 
coefficient of a = 250 ppm/cm while the LMA data indicates a to be in the order of 
;< 5 ppm/cm. Apart from that the two curves show a very similar characteristic. At 
low intensities no non-linear contribution to the free-carrier absorption takes place. 
As outlined in (SJ , from the residual doping of the crystal an absorption coefficient of 
awl ppm/cm can be predicted, excluding this effect as a possible explanation for the 
observed high absorption values. The expected contribution is shown in trace (b) of 
Fig. [3] (solid line). 

A possible explanation for the difference between our observations and the LMA 
data [9] is absorption in the surfaces of the prototype test mass. The crystal surface 
was super-polished to a specified residual microroughness of 2 — 4 A. It is known that 
the polishing procedure can generate thin layers of amorphous silicon (a-Si) on the 



crystal surface 17 . At the wavelength of 1550 nm the absorption of amorphous silicon 
is many orders of magnitude larger than the absorption of crystalline silicon. Although 
the absorption coefficient depends on the thermal treatment of the material, so that 
in contrast to the crystalline material no unique value can be given, absorption values 



of 100 /cm up to 2000 /cm at 1550 nm (0.8 eV) have been reported |18p0] . It is worth 
noting that very few information is available on pure amorphous silicon, since the 
semiconductor industry merely concentrates on hydroginated material (a-Si:H) where 
a vast range of literature is available (see e. g. (2l]). In a different experiment, we 
measured the absorption in a HR stack of Si/Si02 dielectric coatings to a value of 
1000 ppm - 1500 ppm, depending on the polarization of the probe beam. The Si-layers 
of this coating also (mainly) consist of amorphous material. Since in a HR coating 
the penetration depth of the laser beam typically only comprises the first few layers, 
the absorption can be estimated to a value of several tens/cm. This value is even 
smaller than expected for amorphous silicon. However, since thermal treatment and 
the procedure used to apply the Si film strongly infiuence the absorption value (indeed 
a recrystallization can be achieved by tempering) no reliable expectation value can be 



derived from literature. A publication is currently in preparation 22 . For the sake of 
completeness it shall be mentioned that after polishing a layer of native oxide forms 
on the surface and reaches a thickness of 1-2 nm after several weeks-months, strongly 
depending on the environmental conditions [23| . This layer, however, mainly consists 
of Si02 and thus in a first order estimate is not expected to significantly contribute 
to the absorption. 

In view of the above argumentation, for low intensities as planned to be employed 
in ET our measurement results can be interpreted as follows: The total absorption 
consists of 

(i) a bulk absorption ttbuik of a few ppm/cm that is negligible when compared to 

(ii) an offset of about ckoff ~ 250 ppm/cm due to absorption in the crystal surfaces. 
Further measurements for a more exact determination are proposed in section [4] 

Considering the fact that during one cavity round-trip four surface transmissions 
take place (see Fig[I](b)), the absorption per surface can thus be estimated to 
250 ppm/cm x 13 cm/4 « 800 ppm. Except for this offset our data shows a very 
good agreement with the LMA measurements. For comparison. Fig. [4] shows our 
results along with the LMA data to which a constant offset of aoff = 250 ppm/cm 
was added. The very similar characteristics support the assumption of a constant 
offset that cannot be explained by impurities resulting from residual doping. It strikes 
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Figure 4: (a) The data shown m Fig. p^replotted for a comparison to trace (b) taken 
from [9^ that shows the bulk absorption obtained from beam deflection measurements 
on a silicon crystal with a specific resistivity of ~ 30 kf2cm. According to the discussion 
given in the text, an offset of ttofi — 250ppm/cm was added to the data from ref. lo]. 

out that in the LMA measurements the increase of the absorption sets in at lower 
intensities than in our case. As discussed above, a possible explanation for this is the 
higher specific resistivity of the material analyzed by the LMA group that determines 
the lifetime of the generated free carriers. 



4. Conclusions and outlook 

Summarizing, our measurements confirmed the non-linear dependence of the 
absorption in crystalline silicon on the optical intensity reported by Degallaix et al. [9] . 
In addition to the bulk crystal absorption we have observed an intensity-independent 
offset that possibly arises from absorption in a-Si surfaces created during the polishing 
procedure. In a separate measurement 11 it was shown that the absorption in the 



dielectric Si02/Ta205 coatings is orders of magnitude smaller and therefore negligible 
for this consideration. Our results imply an absorption loss of w 1600 ppm for a single 
transmission through a complete test mass, e. g. through the input test masses of the 
ET arm cavities. 

While ET will use large beam diameters and thus rather low optical intensities fi] 
so that the non-linear contributions to the absorption will be negligible, the surface 
absorption might become an issue. The laser power absorbed in the test masses will 
have to be extracted from the system via the suspensions. Because of the small fiber 
diameters only a very limited amount of thermal energy can be extracted. Since 
the total amount of generated heat determines the lower limit for the test mass 
temperature that can be achieved, the high surface absorption will possibly restrict 
the suppression of thermal noise in ET. Therefore, a series of further characterization 
measurements is currently in preparation. In a first step, a set of monolithic cavities 
with different lengths will be prepared and characterized. These results will allow to 
precisely separate the bulk absorption from the surface absorption, since the latter will 
only lead to a constant offset independent from the cavity length. Thereupon crystals 
with a different specific resistivity ranging from 1.5 kf2cm to 70 kJ7cm will be compared 
to analyze the influence of the residual doping on the optical absorption. The goal of 
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these measurements will be to analyze whether crystals produced with the Czochralski 
technique are usable in ET in terms of the optical absorption. Today Czochralski- 
grown crystals that fulfill the ET requirements (450 mm diameter) can already be 
produced. In contrast it is not yet clear whether the considerably purer float zone 
grown crystals will be available in the same dimensions within the next decade due to 
technical difficulties intrinsic to the manufacturing process. Thus, a confirmation of 
the fact that Czochralski-grown crystals fulfill the ET absorption requirements would 
significantly relax the requirements on the ET test masses in terms of availability and 
feasibility. Finally, the discussed measurements will be repeated at the ET design 
temperature of 10 K to characterize the influence of cryogenic temperatures on the 
free-carrier absorption. 
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